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We develop an instanton technique for calculations of correlation functions characterizing statisti- 
cal behavior of the elastic string in disordered media and apply the proposed approach to correlations 
of string free energies corresponding to different low-lying metastable positions. We find high-energy 
tails of correlation functions for the case of long-range disorder (the disorder correlation length well 
exceeds the characteristic distance between the sequential string positions) and short-range disorder 
with the correlation length much smaller then the characteristic string displacements. The former 
case refers to energy distributions and correlations on the distances below the Larkin correlation 
length, while the latter describes correlations on the large spatial scales relevant for the creep dy- 
namics. 
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An elastic string in two-dimensional random potential 
is an archetypal problem of statistical physics, with ap- 
plications to a widest variety of systems and phenomena. 
The examples include vortices in superconductors 111 , dis- 
locations j|] , domain walls Q , and non-Hermitian quan- 
tum mechanics [Q. The string that energetically disfa- 
vors overhangs (like vortices) maps onto directed polymer 
M. In the random potential the string bends to adjust 
itself to the pinning potential relief. Directed polymer 
approximation corresponds to weak pinning, where the 
string shape is smooth and coincides with the path of 
the classical particle traveling through the corresponding 
rugged energy landscape. The Hamiltonian of the string 
has the form flrJ 



H 



dx 



+ V(x,0 







1 


dx 



< 1. 



(1) 



Here x is the coordinate along the preferential direction 
(for small displacements it also measures the length of the 
segment), C( x ) is transverse displacement, k is the elastic 
constant, and the random potential V(x, £) is routinely 
assumed to be Gaussian, with zero average (V) = and 
the correlation function 

(V(x, QV(x', C')> = 0S(x - x')K(C -O, / K(y)dy = 1. 



The amplitude /3 depends on the type of disorder. For 
instance, for a vortex pinned in a plane by point defects, 
(3 w ^niV 2 , with £oi n i, and V being the core radius, 
concentration of defects, and the energy of vortex-defect 
interaction, respectively. 

Much of the early effort was concentrated on the fluc- 
tuations of the free energy of the string. Let the left end 
of the string be fixed at Q = 0, xl — * — oo. The quantity 
of interest is then the free energy of the string e{x, Q as 
the function of the position of its right end (xt, + x, £). 



Knowledge of the mean free energy (e) is important for 
calculation of dynamical quantities, such as pinning en- 
ergy or drift velocity. Various tools, which, among others, 
include numerical analysis |@] , Bethe Ansatz solution (Jt| , 
and power counting based on the effective action , yield 
the result for the disorder-averaged free energy (e) oc x 1 / 3 
and the average displacement of the string (£) oc x 2 / 3 . 

The dynamic response, however, requires more de- 
tailed statistical description, since in such a complex 
system a non-self-averaging behavior can be expected. 
Many of the results have been already made available in 
the course of 15-year-long research, see Ref. |B for an 
excellent review. For recent developments, we cite cal- 
culations of the free energy cumulants || and tails of its 
distribution function |Io| ] . 

So far, attention has been mainly drawn to statistical 
properties of the free energy for a fixed position of the 
end, e(x, £). Of interest are also correlation functions de- 
scribing different positions of the string important for its 
slow dynamics (noise, velocity correlations). Conceptu- 
ally, they also test non-Gaussian shape of distributions. 
In this Letter, we make the first step in this direction and 
evaluate the distribution function, 



P(u) = (S [u-e(0,Ci) + e(0,C2)]). 



(2) 



which describes correlation between the sequential ener- 
gies of the string as it is displaced transversely to the pref- 
erential direction. In what follows, we assume £i > £2- 

Model. The starting point of our theory is the equa- 
tion for the free energy e which can be derived by the 
transfer matrix method from Eq. (Fy) and reads 



de 
d.r 



1 

2k 



de 



T d 2 e 



V(x,0, 



(3) 



where T is temperature. Eq. (|^) is easily recognized as 
Kardar-Parisi-Zhang (KPZ) equation in 1 + 1 dimensions 



1 



M. To this end, x and C can be identified as time and 
coordinate, respectively. 

We focus on the behavior of the correlation function (||) 
at large u. In this case, P(u) is the probability that the 
energy changes considerably upon small displacement of 
the string. This probability is expected to be exponen- 
tially small, and an appropriate technique to calculate 
exponential tails is the instanton approach, allowing us 
to avoid difficulties associated with the replica method. 
The Laplace transform of the function (||) can be pre- 
sented as the functional integral S, 



n(A) 



3 A[e(0,Ci)-e(0,C2)] 



DeDipe 



(4) 



with the effective action 
S[e,i/}] = i I dxd(ip(x,() 



de 1 ( de 



2k d( 2 



-p j dxd(dc^x, c)K(( - ?Mx, a 

+ A[e(0,Ci)-e(0,C 2 )]. (5) 
To arrive to Eq. (@), we have used the identity 



Ded 



®L JLf— V-— — -V 



= const, 



which is a consequence of causality, and subsequently 
performed averaging over the random potential V. 

The field theory (|J) is essentially two-dimensional. To 
move further, we extend the technique developed by Gu- 
rarie and Migdal [[flj (GM), who studied velocity cor- 
relations in the Burgers equation. This is, in order to 
describe the tails of the correlation function (||) we have 
to find the instanton (saddle-point) trajectory and calcu- 
late the action $,•„ at this trajectory. The result for 11(A) 
with the exponential accuracy reads 



11(A) =exp [S in (X) - S in (0)] 



(6) 



Prefactors can be obtained systematically by expanding 
the action around the instanton path, but this goes be- 
yond the scope of this Letter. 

The saddle-point equations for the effective action 
S[e, ip] have the form 



de_ 
dx 



1 

2k 



2k d( 2 ~ 
T d 2 ip 



dip l^_f,2l\ 

~dx~ + Kdc \ <KJ + "2k de, 2 



-20 j dCK((-c)i>(xx'y, 

= -i\6(x) [6(( - d) - 5(( - ( 2 ) 



(7) 
(8) 



Long-range disorder. We assume first that the func- 
tion K is only slightly changed on the scale of Ci — C2- 
It can be then expanded, K(y) = ko — k\y 2 j2. In the 
subsequent analysis, we follow GM. 

Note that Eq. (|J) can only have non-zero solutions for 
x < ("diffusion in reverse time"). The field ip is dis- 
continuous at x = 0; integrating Eq. (||) between —0 and 
+0, we obtain the boundary condition, 



V(s = -o,C) = *M*(C-Ci)-*(C-C2 



(9) 



To solve the saddle-point equations, we disregard the 
term with d 2 ip/d( 2 in Eq. (|J) for a while. Then, due to 
Eq. (||), the function ip always remains the sum of two 
delta- functions, 

i{,(x, C) = »/i(x) [S(( - Pl (x)) - S(C - P2(x))} , (10) 

where the tilde means that this solves the "incomplete" 
equation. Eq. @ implies p,(0) = A, pi,2(0) = Ci.2- From 
Eq. (Q) we then find e(x, £) = d(x) + 6(x)£, with 

d' + 1 2 /2k = -Phuipj - pi); V = 2f3k 1 p( Pl - P2 ). 

Substituting this back into Eq. (|^), we find p{x) = X 
and p\ (x) — P2 (x) = Ci — C2 • These solutions and, conse- 



quently, the resulting form of the distribution function, 
are very different from those for Burgers equation fll"|| . 

Now we include the term with d 2 ip/d( 2 into consider- 
ation. We look for a solution of Eq. (]?]) in the form of 
a linear function of £, e = a(x) + 6(x)£ (the consistency 
is checked afterwords) . The linear equation (^) is easily 
solved and yields 



iXK 1 ' 2 

(2 7 rT|x|) 1 /2 



"(C-P2) 

e 2Ta: 



(11) 



with p\ = p' 2 = b(x) j k. Substituting this into rhs of 
Eq. (0), we find that the latter is temperature indepen- 
dent, i.e. it remains the same as that without the dif- 
fusion term. Thus, the linear Ansatz for e is consistent. 
Also, we find b' — 2(3k\X{pi — p%)- The solution is then 

Pl,2(x) = Cl,2 + b Q X. 

Next, we calculate the action S along the instanton 
trajectory. As seen from the saddle-point equations, the 
first term in Eq. (||) is just doubled the second one. Mul- 
tiplying Eq. (^J) by e, integrating it over x (from —00 to 
—0) and over £, and comparing with Eq. (0), we find 
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A[.r((U ,!---((>.,, I] =-2.1 j vKv~J dxdQv 



Since e depends on C, linearly, and J ipdQ = 0, the last 
term in the right-hand side vanishes. The instanton ac- 
tion acquires the form 

S m = -(3 ( dx f d(dC^(x,OK((-C)i>(x,('). (12) 



Substituting the instanton solution (|11| ) into Eq. (|12|) , we 
find that the term proportional to fco i n the effective ac- 
tion vanishes, while the contribution with k\ diverges in 
the limit of large negative x. This is because our consid- 
eration is limited to distances shorter than the correlation 
length of the random potential £ ~ (ko/ki) 1 ^ 2 , so that 
we could replace K(y) by its expansion. This cut-off ef- 
fectively defines the constant bo- Replacing integrals over 
C in Eq. ( |l2| ) in infinite limits by integrals from — £/2 to 
£/2 and calculating them, we arrive at the cut-off for x, 



-k£ 2 /T, 

-«//3fc!ACo) 1/2 , 



T>T C 

T <Tn 



(13) 
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FIG. 1. Distribution function P(u) for long-ranged corre- 
lated disorder, Eq. ([li]) . 

with Co = Ci - C2 > and T c = (/3fciKC £A) 1/2 . The in- 
stanton action reads Si n = /3&iCo A 2 |x c |. Note that, since 
our cut-off procedure is somewhat arbitrary, we have re- 
moved all numerical factors from the action. Laplace- 
transforming Eq.(^|), we arrive at the expression 



P(u) 



(14) 



which is illustrated on Fig. 1, with c ~ c' ~ 1. 

The result ( pT| ) was obtained by instanton approach, 
which means it has to be exponentially small. In partic- 
ular, this requires u > 0. We thus get the conditions u 3> 
CoiPkint;) 1 / 3 for the elasticity-controlled regime (P(u) oc 
exp(-u 3 )) and Co(/3fci< 2 /T) 1/2 < u < Co?7£ for the 
temperature-controlled regime {P(u) oc exp(— u 2 )). The 
latter result only holds for temperatures higher than 
Td = £ 4 / 3 (/3fciK) 1//3 , which does not depend on Co and 
has a meaning of depinning temperature — the typical 
pinning energy on the Larkin correlation length. Above 
the depinning temperature correlation functions decay 
faster than exponentially Q. 

Short-ranged disorder. We take now K(y) — 6(y). 
At x = 0, the field ip is a set of a delta-peak (located at 
C = Ci) an d a delta-dip (£ = (2), see Eq. (||). As we 
trace the evolution in reverse time, these features smear 
and move. For short times, before (if ever) they intersect 
or smear so much that they become indistinguishable, ip 
remains a peak-dip function, centered at £ = pi(x) and 
£ = pi{x), and sharply vanishing away from these points. 

Let us calculate the energy created by a peak-shaped 
field if>, located at £ = p, far from this point. At large 
distances from the peak, the solution must have a scal- 
ing form. Eq. (Q) with the zero right-hand side allows 
for only one scaling solution e = /((£ — p)/x 1 ^ 2 ) vanish- 
ing at x — > -co, which at large scales, where the term 
with the second derivative can be disregarded, becomes 
e = k{( - p) 2 /2x. For |C - p\ < {Tx/k) 1 ' 2 the diffusion 
term is important, and this simple form of the scaling 
solution does not apply. 

Our solution has an obvious drawback: It diverges as 
x goes to zero. To amend this, we write 



<C-p) 2 

2(x - x ) ' 



x > 0. 



(15) 



The regularization constant xq plays an important role 
and is found from the following considerations. For 
x = 0, the energy given by Eq. ( |T5| ) is e — — K C 2 / x o- I n 
the scaling regime, the only one relevant length scale is 
Co = Ci — C2 ; and the only relevant energy scale is tem- 
perature. Therefore, T ~ k(q/xq, whence xq ~ k£q/T. 
The numerical factor remains undetermined, but can be 
deduced from numerical solutions of Eqs. ([?]), (§h. 

Now we return to Eq. (||) and solve it at small |ac|, 
when the peak and the dip are well separated. We write 
ip(x,Q in the form ip pe ak + ipdip- The peak (located 
near £ = pi(x)) sees then the energy produced by it- 
self and the energy produced by the dip. The energy of 
an isolated peak (|l5|) grows as C 2 far from the peak, and 
hence at the point p\ the energy produced by the dip is 
much greater than the peak contribution. If, in addition, 
P2 — Pi ^ (T(xq — x) /k) 1 / 2 , we can use Eq. ( ffa ) to eval- 
uate the coefficient de/d( in Eq. (^). We get then the 
equation 
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dlppeak 4>peak Pi ~ P2 (typeak T 9 2 (f> 

dx x — xq x — xq d( 2k d( 2 



0, (16) 



with the boundary condition tppeak{x — 0) = iA<5(C — Ci)- 
Writing a similar equation for the dip (C « pq) and 
solving both equations, we obtain 



Ippeak 

P'l 
P2 



= ± 



± 



iXxon 1 / 2 



1 



(2wT\x\) 1 / 2 Xq _ x 

Pi - P2 
x — xq 



exp 



«(C-pi, 2 ) : 



2Tb 



(17) 



As follows from Eq. ( |l7[ ) , positions of the peak and the 
dip are pulled apart. The distance between them grows 
as (x — xq) 2 ; at the same time, they smear as |cc| 1//2 ; thus, 
they never overlap and can be considered as well sepa- 



rated at any x. Using Eq. (17) to calculate the instanton 
action, we get 



a(3X 



2 K Co 
T ' 



Co - Ci - c 2 



1. 



(18) 



which translates into the expression for the distribution 
function 



P{u) = exp 



aTu 2 



(19) 



with a numerical constant a of order one. The expression 
( |l9| ) is valid for temperatures much higher than /3k£o /u 2 
and is an analog of the upper line of Eq. @. We conjec- 
ture that the long-u tail of the distribution function (the 
analog of the lower line of Eq. (|l4|)) also exists, though 
we were not able to obtain it with this method of solving 
Eqs. §, ®. 

Discussion and conclusions. Eqs. (jlj) and jl9| ) 

constitute the central result of this Letter. The function 
P(u) has a meaning of the probability that the differ- 
ence of energies of a string of the same length and the 
transverse displacements Ci and C21 equals u. An instan- 
ton solution corresponds to the situation when at Ci is a 
minimum of the potential energy while at C2 is the neigh- 
boring maximum. P[u], thus, measures the height of the 
barrier of the free energy relief. For large u to find a 
high barrier is quite improbable, consequently the result 
is exponentially small. 

If the string has a finite transverse size £0; the correla- 
tion radius £ of the disordered potential is of order of £o- 
Eq. (|l4|) describes thus the case when the string is dis- 
placed at a distance small compared with its transverse 
dimension, and Eq. ( |il^ ) applies in the opposite regime 
of long distance. 

These results have an important conceptual value, 
since they describe correlation functions of energy at dif- 
ferent positions for the directed polymer problem. We 
argue now that they also provide certain predictions for 
experimentally observable dynamical properties of the 



string. We suggest two types of experiment, which cer- 
tainly do not exhaust all the opportunities. 

(i) Creep of domain walls. If a domain wall moves 
from one position to another one, the distance between 
the two positions can be measured. The quantity Coj 
which is the distance between a position and a maximum 
separating the two, is not known, but can be determined 
if the wall travels back and forth. The tunneling rate 
is determined by the height of the barrier u. A similar 
information can be extracted from the random telegraph 
noise resulting from the motion of the string between two 
positions. The time spent in each position is determined 
by the barrier to be overcome seen from this position. 
These considerations do not take into account a possibil- 
ity of macroscopic quantum tunneling, which has to be 
considered separately. 

(ii) Distribution of pinning energies. The string is (lo- 
cally) depinncd if it moves from one free energy mini- 
mum to the next one. The pinning energy corresponds 
precisely to the quantity u above. However, globally one 
has to average over the distribution of the displacements 
Co, i.e. to know the distribution of distances between 
adjacent minima and maxima of the random potential. 
Measurements of the distributions of the pinning energy 
compared with our results can provide an information 
about this distribution, which, to our knowledge, has not 
been previously discussed. 

In conclusion, we have developed an instanton ap- 
proach to calculations of various correlation functions 
describing statistical behavior of the elastic string in the 
two-dimensional disordered potential. We applied our 
technique to the investigation of correlations of free ener- 
gies corresponding to different low-lying metastable posi- 
tions of the string. We have found the asymptotic behav- 
ior of such energy-energy correlations for the moderate 
spatial scales (within Larkin correlation length) and the 
large scales, exceeding Larkin length. The latter situa- 
tion corresponds to the conditions of the creep dynam- 
ics. We have discussed applications of our results to the 
dynamic response and noise in various two-dimensional 
systems such as domain walls, vortices, and dislocations 
in thin films. 
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